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TECHNICAL NOTE NO. 1808 


METHOD FOE EVALUATING FROM SHADOW OE SCHLIEEEN 
PHOTOGRAPHS THE PRESSURE DRAG IN TWO-DIMENSIONAL OE 
AXIALLY SYMMETRICAL FLOW PHENOMENA WITH 
DETACHED SHOCK 
By Antonio Ferri 

SUMMARY 


A method has "been developed for evaluating from shadow or schlieren 
photographs the pressure drag of axially symmetrical todies at zero 
angle of attack or of two-dimensional todies producing detached shock. 

The method consists in the determination of the flow properties along 
a characteristic line and the application of the momentum theorem to a 
stream tute around the tody. The method can te applied to the determina- 
tion of external drag of supersonic inlets with subsonic flow at the 
entrance of the Inlet. 


INTRODUCTION 


The system of characteristics permits the pressure drag and the 
shape of the shock for axially symmetrical and two-dimensional flow 
phenomena having everywhere supersonic speed to te determined directly 
from the boundary conditions. In this case if viscous effects are 
neglected, the value of the pressure drag can te determined analytically, 
and no experimental measurement is necessary. The determination of the 
shape of the shock can te useful in making a comparison between experi- 
mental and theoretical results and permits an evaluation of the displace- 
ment thickness of the wake and of the boundary layer along the tody. 

The characteristic system does not permit, however, direct determination 
of the pressure drag when the shock is detached from the tody, because it 
cannot te applied to the subsonic part of the flow. 

The determination of the pressure drag of todies producing 
detached shocks is important in many problems; for example, the todies 
of revolution used in practical applications often have round noses and, 
in supersonic inlets, conditions exist in which the flow at the entrance 
of the inlet is subsonic. The experimental measurement of the pressure 
drag for these cases usually requires complicated experimental equipment; 
therefore, the determination of drag from shadow or schlieren photographs 
of the phenomena can te very useful for practical applications. A 
method is presented herein that permits determination of the pressure 
drag directly from shadow or schlieren photographs. This system can te 
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especially useful in all cases in which, the experimental technique uses 
the principle of moving models (shooting ranges, whirling arms, etc.) 
and in experiments of supersonic inlets in which also the internal flow 
produces forces and in which it is, therefore, difficult to measure the 
external drag independently of the internal drag. The method gives the 
value of the pressure drag due to the subsonic and supersonic regions of 
the flow and permits also the determination of pressure and velocity 
distribution of the supersonic region, from which the shape and thickness 
of the wake can be obtained. 


SYMBOLS 


X, 7 

Cartesian coordinates 

n 

normal, to streamline 

V 

velocity 

Y l 

limiting velocity 

W 

velocity in terms of ; the limiting velocity (V/VT ) 

M 

Mach number 

P 

density 

P 

pressure 

s 

entropy 

m A , mg 

mass flow 


sin p tan p sin cp „ . _ 

l = for axially symmetrical phenomena; 1 = 0 for 

cos(<p + p) 

two-dimensional phenomena 


m = 


sin p tan p sin <p 


cos(cp — p) 

two-dimensional phenomena 


for axially symmetrical phenomena; m = 0 for 


Inclination of the velocity vector with respect to x-axis 


€ 


inclination of the shock with respect to upstream velocity 


Mach angle 



\ 

* 




P 
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c T specific heat at constant -volume 

7 ratio of specific heats 

R gas constant 

Subscripts-: 

0 refers to stagnation free-stream. conditions 

1 refers to free-stream properties 

A points of first family 

B points of second fam il y 

C quantities at the points calculated from. A and B 

P quantities at any point P "behind the shock 

THE METHOD FOR DETERMINING THE PRESSURE DRAG 


In any supersonic phenomenon in which a detached shock exists, the 
flow is subsonic only in a small zone of the space and the speed outside 
of this zone, which is usually limited to the vicinity of the nose of 
the body, becomes again supersonic. The body considered has finite 
dimensions. In' the supersonic part of the flow the characteristic theory 
can be applied to determ i ne the pressure drag and the pressure distribu- 
tions. The line BN of figure 1 is the sonic line that divides the 
subsonic from the supersonic region, and from any point G of the shock 
in the supersonic region the characteristic line (3) of the second 
family that meets the body at some point D (see, for example, reference l) 
can be determ i ned. 

If the flow properties are known along the line (3), the component 
of the momentum of the stream along (3D in the direction of the free- 
stream velocity and the resultant of the pressure on the surface GD in 
the same direction can be evaluated. The resultant of the pressure 
along the surface GD in the direction of the x-axis can be obtained 
from the momentum law. Indeed, the difference of momentum between the 
surfaces EE’ and GD must be equal to the resultant of pressures at 
the surfaces E’G, EE*, E0, 0D, and GD. The pressure along E*G and E0 
does not produce any component of force in the direction of the x-axis. 

The resultant of pressure at the surfaces EE* and (3D is known (along IE* 
the pressure is equal to the free-stream pressure); therefore, the 
resultant of the pressure in the direction of the x-axis along 0D can 
be evaluated from the difference of momentum in the direction of the 
x-axis between EE* and (3D. "When supersonic Inlets having subsonic 
entrance velocity are considered (fig. 2), the streamline E0 is not 
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parallel to the stream direction and, therefore, with this system the 
pressure drag along this line is also determined. 

In these considerations the viscous effects have teen neglected, tut 
the viscous effects are small 1 in the front part of the tody and therefore 
can he considered independently. For the tack part of the tody the 
considerations made in the "Introduction" on the possibility of the 
determination of the displacement thickness of the boundary layer and of 
the wake are still valid. It should be noted that, by consideration of 
the change in momentum in the y direction, the lift of two-dimensional 
todies at an angle of attack can be determined. From the preceding con- 
siderations the determination of the pressure drag along EQD (fig. l) iB 
reduced to the determination of the flow properties along GD. The flow 
properties along (3) can be determined from the shape of the shock in 
the following way: 

If EGL is the shape of the shock determined from experimental 
measurements (shadow or schlieren photographs), the inclination of the 
shock along GEL (fig. l) can be measured. In order to obtain the inclina- 
tion of the shock with same precision the coordinates of the shock can 
be measured and an analy tical expression can be determined for the curve. 

From the inclination of the shock 6, the magnitude W and the 
direction cp of velocity and the variation of entropy Ab across the 
shock can be calculated at any point P of the shock, for example, from 
the equations: 



7 + 1 % 2 

2 M L 2 sin 2 e - 1 



tan e 


( 1 ) 


tan €p 
tan( e - <p) p 


7-1 
7 + 1 


1 — Wp 2 cos 2 (e — qp) 


Wp 2 sin 2 (e 


9) 


( 2 ) 


Pp tan 6p 
P-L tan( e - cp) p 


( 3 ) 
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27 

7 + 1 



sin^e 


P 





ASp = Cy logg 


Pp/V 

p l\ p p, 


(5) 


In the equations, the subscript P indicates quantities at the 
point P behind the shock and the subscript 1 indicates free— stream, 
conditions. Prom two points of the shock A and B (fig. l), the 
flow properties at the point C can be determined by means of the charac- 
teristic system for rotational flow as is shown in the following 
discussion. 

From A the tangent to the characteristic line of the first family 
and from B the tangent to the characteristic line of the second family 
can be drawn (reference l) 



tan(cp + 



> 



= tan(9 — p) 


B 


(6) 


where 


Bin^p 


7 - if 1 

2 W 2 


1 


The intersection of the two characteristic lines determines the 
point C in the first approximation. The properties at C in the first 
approximation can be determined from the following expressions (reference l) 
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%~^C 

¥ 


— tan 




X A ~ *C 
7 


1 


de _1_ sln3p / 
dn 7R cos(p + qp) ' ^ 


Xq) = 0 (7) 


W B ~ ¥ C 
¥ 


+ tan 





ds 1_ sln3p / 
da JR cos(<p — p)^® 


*0) = 0 ( 8 ) 


where, for two-dimensional phenomena, 2 and m are zero whereas , for 
axially symmetrical phenomena. 


sin p tan (3 sin cp 
cos(cp + p) 


(9) 


sin p tan p sin qp 
m = — 

cos(qp — p) 


( 10 ) 


la the first approximation the quantities with no index in equation ( 7 ) 
can he assumed equal to the come speeding quantities at the point A 
rnifl in equation ( 8 ) to the quantities at the point B, and for ds/dn the 
following expression can he used: 


ds 

da 


^A 

sin p A / \ sda Ps 

' -A Z 0'cos(p + cp) A + I ^'cos(qp — p)p 


( 11 ) 


Tn this approximation the entropy gradient is assumed to he the same 
in the zone between B and C as in the zone between A and B. In 
the process of calculations it is possible that the two Known points 
fall on the same streamline. In this case, the value of ds/dn for the 
first approximation can he assumed equal to the value of ds/dn at one 
of the two points. The value of ds/dn can he obtained from the differ- 
ence in entropy at this point and at a nearby knova point on the same 
characteristic line. 

After the properties at the point C in the first approximation are 
obtained, a second approximation can be obtained in the following way: The 

value of entropy at the point C is determined from the entropy of the 
point on the shock that is in the same streamline as . the point C. This 
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point is determined from mass— flow considerations „ The flow at C is 
then determined again by averaging the values of the known terms in 
equations ( 6), ( 7 ), and (8) between the values at A and C at B 

and C, respectively. 

For axially symmetrical phenomena the mass flow contained in the 
stream tube bounded by the two surfaces of revolution having as generatrix 
the streamline that passes at A and C (m^) and at B and c (mg) can 
be determined from figures 1 and 3 : 


m. 


J ?c 




A sin B 

A 2*ypY ; — - — r dy 

sin(cp + p) 


mg 


< 


B sin p 

-2irypV — dy 

sin(cp - p) 


( 12 ) 


The density at any point P can be expressed as a function of the 
free-stream stagnation density and of the local velocity and entropy in 
the form 


Pp - 



(13) 


and, therefore, can be calculated at A, B, and C. The velocity at any 
point along AC and BC can be expressed in the form: 


Y = 


1 + 


(j ~ 7 A. 




Therefore, 


t = t a§ + *lO-Ja)] 
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Similarly, 


sin Pg 

“B = rtP B Y B sinfo' - p) B ( ' 7 C “ 



+ 


2 y c + y B 

3 




+ B r 



■where 


- 



1 

?C ~ 7 3 




Bg = 




1 

~ fb 




C 2 ~ 


sln Pc sin((p - p) B _ 
sin pg sin(cp - p) c 


7C ~ 7B 


If a two-dimensional phenomenon is considered. 


m A = -PA T A 


sin Pa wn J h 1 /" 

ia(cp + p). (7 C ” 7 k ) + 2 C A 1 + B 1 + C i 


( 7 c ~ 7 a 1 


\2 




sin ^ 
sin(qp — P) 


B 


(y c " 7 b) + ^ c P ~ ~ (h + B 2 + c 2^ 


(16) 

(IT) 

) (18) 

l (19) 
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When the mass flow is determined, the points C* and C** on the 
shock along the same streamline that passes at C can he dete rmin ed for 
axially symmetrical phenomena from the expressions 


m A = p i y i(7a 2 “ r c 2 > 


“B = Pl Y l(y-B 2 - To 2 )* 
and, for two-dimensional phenomena, 

®A = P^Cta - 7 q ) 
“B = “ ^c) 


In the first approximation, C 1 will he slightly different from C* *; 
therefore, a point C 1 ** will hare to he chosen between the two points. 

The variation of entropy at C*** can he determined from equation (5), 
and the entropy at C ,,f behind the shock iB the same as the entropy at C„ 
Wow at the point C, the quantities As, <p, p, and W are known in the 
first approximation; therefore, a second approximation can he determined 
by ass uming for the quantities without index in equation ( 7 ) the corre- 
sponding average values between A and C, and in equation (8) the average 
values between B and C. In equation ( 7 ), 


ds _ 
dn “ 


( S A “ So)Jcos((p + p) A + cos(9 - pjJ 

( z a - ^cX 8111 Pc + sl11 Pa) 


and in equation (8), 


( S B “ s c) cos (<P - P)b + cos(cp - p)qj 
(*B ~ *cX Bla P C + sijl Pb) 


ds 

dn 
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A second approximation can "be obtained also far the points C* 
and C** and, therefore, for the value of As at C. If necessary, 
a higher approximation can he calculated for As. Proceeding in a 
similar way permits all, the characteristic net CDL to he determined 
and, therefore, the flow properties along CD can he obtained. For 
inl ets the calculation also gives the value of mass flow that goes 
inside the inlet, because it gives the position of the point , E that 
limits the stream tube that goes inside the Inlet. 

A check of the precision of the calculations can be obtained from 
the comparison of the mass flow contained in the stream tube EE* with 
the mass flow across CD. In addition, the entropy at D must be equal 
to the entropy at 0 (fig. l). 

If the tests are performed on a fixed model, a pressure measurement 
at a point in the zone DF can also be used for control of the precision 
of the calculations. This control is useful when inlets are considered 
because the stream tube E**E* (fig. 2) in this case is not known. 'When 
the flow field along CD has been determined, the pressure drag along EOD 
can be determined from the momentum equation, and from the shape of the 
shock GL and from the characteristic line GD, the pressure along DF 
and the shape of the streamline DT can be determined. The streamline DT 
gives the displacement thickness of the boundary layer and of the wake. 


conchjsions 


A method for determining pressure drag in axially symmetrical or 
two-dimensional flow phenomena from shadow or' schlieren photographs has 
been developed. The method uses the characteristic system for rotational 
flow and can be applied when detached shocks are produced by the body. 

The method can be applied to the determination of mass flow and 
external drag of supersonic inlets with subsonic flow at the entrance of 
the inlet . 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., December 22, 19W3 


REFERENCE 


1. Ferri, Antonio: Application of the Method of Characteristics to 

Supersonic Rotational Flow. MCA Rep. No. 84l, 19^6 . 






